Meissner effect cannot be explained classically 
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The Meissner effect is an important characteristic of superconductivity and is critical to distinguishing super- 
conductivity from simply the absence of electrical resistance (perfect conductivity). In a recent paper published 
in American Journal of Physics, Essen and Fiolhais claimed that the Meissner effect is explained by classical 
physics. [Am. J. Phys. 80 164, (2012).] We claim it cannot be understood by classical mechanics and point 
out that their derivation of the Meissner effect by classical physics is based on an inadequate treatment of the 
magnetic field energy. A correct treatment of the magnetic field energy clarifies the need for quantum mechanics 
to understand the Meissner effect. We stress that Meissner effect is energetically favorable due to the energy 
of condensation of the Cooper pairs. The condensation of electrons into Cooper pairs is best understood as a 
quantum mechanical phenomenon. 
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I. INTRODUCTION 

When we apply a magnetic field to a normal metal, the mag- 
netic field penetrates the metal. If the metal becomes super- 
conducting when the metal is cooled, the magnetic field in 
the bulk of the superconductor may be expelled. The mag- 
netic field is expelled, if the strength of the field is lower 
than temperature-dependent critical field of the superconduc- 
tor This expulsion of the magnetic field is known as the 
Meissner eff'ect.' This phenomenon is distinct from that of 
perfect conductivity. In a classical system with zero resistance 
Lenz's law guarantees that the magnetic flux density cannot 
change in time. Thus if a system becomes a perfect conduc- 
tor below some critical temperature, its magnetic flux density 
will depend on its history; i.e. whether the external field was 
applied before or after cooling. For sufficiently weak fields, 
superconductivity is a well-defined thermodynamic phase in 
which the equilibrium state of the system is independent of 
history. The expulsion of the field contradicts with the Lenz 
law of the classical electrodynamics, and cannot be explained 
without quantum mechanics. The expulsion of the magnetic 
field costs energy which is supplied by quantum mechanical 
condensation energy of the electrons into Cooper pairs. 

It is also noteworthy that superconductors are classified 
into two types. For type I superconductor, there is only one 
temperature-dependent critical field, HdT). If the strength 
of the magnetic field exceeds this value, the superconductiv- 
ity is destroyed, and the magnetic field penetrates the bulk. 
A schematic phase diagram of the type I superconductor is 
shown in Fig. 1. For type II superconductor, there are two 
critical fields, Hci(T) and Hc2{T). When the magnetic field 
satisfies Hc\(T) < H < HdiT), quantized magnetic fluxes 
penetrate the superconductor, forming an Abrikosov lattice.^ 
The magnetic flux is quantized in units of the elementary flux 
quantum hjle, where h is the Planck constant and e is the el- 
ementary charge. This quantization clearly indicates that this 
phenomenon originates from the quantum mechanical effect. 

The magnetic field is expelled from the superconductor via 
a dissipationless supercurrent that flows at the surface of the 
superconductor This supercurrent flows such that the mag- 
netic field created by this current just cancels the external 
magnetic field in the bulk of the superconductor The su- 
percurrent is possible because of a complex order parameter 
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FIG. 1. A schematic phase diagram of a type I superconductor. 
Strength of the critical magnetic field H^.(T) is shown by a thick line 
as a function of temperature T . In the superconducting state below 
the line, applied magnetic field is expelled from the superconductor 
so that in the bulk of the superconductor the magnetic flux density B 
is zero. In the normal state, applied field penetrates the metal. The 
transition temperature to the superconducting state in the absence of 
a magnetic field is shown as T^. 



which describes the macroscopic wave function of the con- 
densed Cooper pairs. The phase of the order parameter must 
be single valued in the bulk of the superconductor This condi- 
tion makes dissipationless current in the superconductor pos- 
sible in spite of the presence of impurities and phonons which 
give finite resistivity in the normal metal. The realization 
of the supercurrent can only be explained quantum mechan- 
ically. The relation between the complex order parameter and 
the supercurrent is supported by much experimental evidence 
including the Josephson effect.^ 

The existence of the supercurrent, however, is not enough to 
explain the Meissner effect. When a system crosses the phase 
boundary shown in Fig. 1 from normal side to superconduct- 
ing side, the supercurrent begins to flow spontaneously, and 
the magnetic field is expelled. In this process the magnetic 
flux in the superconductor changes from finite value to zero. 
According to classical electromagnetism, Lenz's law works to 
resist the change in the magnetic flux. To expel the magnetic 
field from the bulk of the superconductor against this law, it 



is required that the superconducting state in which the mag- 
netic field is expelled has lower free energy than the state in 
which the magnetic field remains in the bulk. In the absence 
of the magnetic field, the superconducting state has lower free 
energy than the normal state provided T < Tc, where T^ is the 
transition temperature to the superconductivity. In the pres- 
ence of the magnetic field, superconducting state that expels 
the magnetic field competes with the normal state in which the 
magnetic field remains. Superconductivity cannot be realized 
when the magnetic field penetrates the bulk, so the supercon- 
ductor needs to expel the magnetic field. However, expelling 
the magnetic field costs energy, so when the applied magnetic 
field becomes strong enough, the energy of the superconduc- 
tor becomes higher than that of the normal state, and the su- 
perconductivity is destroyed. This is the origin of the critical 
field. 

The lowering of the energy in the superconducting state 
comes from the condensation energy of the electrons into 
Cooper pairs. It is described by the BCS theory,^ which is en- 
tirely based on the quantum mechanics. The temperature de- 
pendence of the critical field observed experimentally is quan- 
titatively explained by the theory. Thus, quantum mechanics 
is indispensable for the explanation of the Meissner effect. 

In spite of all this, Essen and Fiolhais^ claimed that the 
Meissner effect can be explained by classical mechanics.* 
From our argument above, it is evident that their claim is in- 
correct. It is however instructive to examine their claim and 
understand the flaw of their argument. In the next section, 
we introduce their argument. In Sec. Ill, we see how several 
strange conclusions are derived if we accept their argument. 
In Sec. IV we clarify the flaw of their argument and find that 
they treat energy of the magnetic field incorrectly. We need 
only classical electromagnetism to understand the flaw in their 
argument. No understanding of quantum mechanics or super- 
conductivity is necessary to understand our arguments in this 
section. 



II. ARGUMENT BY ESSEN AND FIOLHAIS 

Essen and Fiolhais try to persuade us that Meissner effect is 
a property of a perfect conductor So let us forget the fact that 
the supercurrent can only be explained quantum mechanically, 
and try to examine whether or not a perfect conductor shows 
the Meissner eff'ect. Their argument is based on a description 
in chapter 1 of a famous textbook by de Gennes,^ where the 
London equation is derived using only classical mechanics. 
The derivation is not appropriate as we will see in Sec. IV, but 
we repeat the argument here. 

For simplicity, the discussion will be Umited to T = 0. It 
is assumed that total energy of superconductor in a magnetic 
field is given by 



E — Es + Eg + Eji, 



(1) 



where Es is the condensation energy of electrons, which 
comes from quantum mechanical effects. The second term. 
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is the energy of the magnetic field applied externally with //q 
being permeability of the vacuum,**, B = jj.oH is the magnetic 
flux density with H being the applied magnetic field and 
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(3) 



is the kinetic energy of the (super) current. In this equation m 
is the mass of an electron, n{r) is the number density of elec- 
trons, and v(r) is average velocity of electrons. This energy 
can be rewritten as follows. The supercurrent density j(r) is 
expressed as 



j(r) = n(r)ev(r) . 



(4) 



It is also related to rotation of the magnetic flux density owing 
to Ampere's law that is valid for a steady current: 

VxB(r)^fioj(r). (5) 

Using these equations the kinetic energy is expressed as 
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(6) 



where A - -\Jm/jJone^ is the London penetration depth. Now 
de Gennes assumes Es to be a constant, and minimizes the 
energy Eg + Ek with respect to the field distribution B(r), and 
obtains the London equation: 



B + A^V X (V X B) = , 



(7) 



(2) 



as a condition to minimize the energy. A solution gives a state 
in which supercurrent flows at the surface layer of the conduc- 
tor, the thickness of the layer being of the order of A,'^ and the 
magnetic field in the bulk is eliminated. This situation is iden- 
tical to the state realized by the Meissner effect. The essence 
of the argument is that the energy of the magnetic field Eg 
is minimized by making B = in the bulk. In this deriva- 
tion quantum mechanics is not used at all. Essen and Fiolhais 
claim based on this argument that the Meissner effect can be 
explained classically. However, this coincidence is superficial 
as we will see later. 



III. STRANGE CONSEQUENCES 

If this argument by de Gennes is correct we encounter var- 
ious strange conclusions. We point out some of them in this 
section. 

(1) In the argument in Sec. II, the condensation energy of 
the superconductor has no role. Thus there is no reason for the 
superconducting state to be destroyed by an arbitrarily strong 
magnetic field. It cannot explain the existence of the critical 
field He above which the superconductivity is destroyed and 
the magnetic field penetrates the bulk. 

Furthermore, unphysical latent heat must exist at the phase 
transition point. Let us consider what happens at the phase 
transition point. We apply a strong enough magnetic field to 
a sample, and cool it to a temperature near absolute zero. The 



magnetic field is stronger than the critical field, so the mag- 
netic field penetrates the sample, and the sample is in the nor- 
mal state. (See Fig. 1.) Then we reduce the magnetic field. 
At H - Hc{T), the sample transitions to a superconductor, 
and the magnetic field is expelled. Let us consider a cylin- 
drical sample with radius R, length /, and volume V - nR^l. 
The magnetic field is applied parallel to the axis of the cylin- 
der We investigate what happens to the energy in Sec. II at the 
transition. Before the transition, there is no current, and the 
magnetic field stays in the sample. The energy of the mag- 
netic field Eg-"^ and the kinetic energy Zs^.'"* in the normal state 
are easily calculated using eq.(2) and eq.(3), respectively. The 
answers are £5*"^ = (l/2/io)fi^V, and E^*"* = 0. 

Next we calculate Eb-^^ and ii^**^ in the superconducting 
state. The supercurrent flows at the side wall of the cylinder 
and cancels the magnetic field in the bulk. The current flows 
in a thin skin layer of the order of A. Since A is typically small, 
of the order of 10"^ m, we can consider that the magnetic field 
is almost completely expelled from the volume V for a sample 
with R^ A. Thus, Eb^'^ = 0. On the other hand, Ei,^'^ is not 
zero in this phase. From the London equation eq.(6), we know 
that the magnetic field decays at the surface as B exp[(r-7?)//l] 
with Q < r < R being the distance from the central axis of 
the cylinder Using this magnetic field, we obtain |V x B| ^ 
{BjA) exp[(r - R)IA]. Then £■['* is calculated by eq.(4) as 
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In the second line we used the fact that the integrand has non 
zero value only at r ^ R. Gathering these results we find that 



energy in the normal state £*"' - E 
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much larger than that in superconducting state /s*^'' = Eg + 

4'' = (l/2;Uo)B^V(/l//?), the ratio being £'"V£*'^ =R/A » 1. 
The difference of these energies and the condensation energy 
£5 that is neglected in Sec. II must be compensated at the tran- 
sition point in the form of a latent heat of the phase transi- 
tion. However, such a latent heat is not observed experimen- 
tally. Actually, the latent heat of a phase transition comes from 
change in the entropy, so it must be zero at T = owing to the 
third law of thermodynamics. The argument of the previous 
section has no experimental support. 

(2) A perfect conductor has no mechanism to dissipate en- 
ergy. Thus there is no process to bring the system to the ther- 
mal equilibrium state in which the free energy takes its min- 
imum value. The state in which Eb = in the bulk is never 
realized, if we start from a state in which Eb > 0. 

(3) In the derivation of the London equation, the fact that 
the metal is in superconducting state is not used at all. There- 
fore, if this derivation is correct, any metal must show the 
Meissner effect. '° 

(4) Let us consider a long cylindrical bulk superconductor 
When the external magnetic field is parallel to the axis of the 



cylinder, supercurrent flows along the surface of the cylinder, 
and the field in the bulk is canceled. This situation is like a 
solenoid used to cancel the magnetic field inside. Now be- 
cause the current flows only within a thin region at the surface 
of the cylinder, nothing will change when a hole is introduced 
along the axis of the cylinder Namely, we will have the same 
situation when we replace the bulk cylinder by a hollow tube; 
Ek stays the same and Eb is independent of whether the space 
is filled with metal or not. However, experimentally the bulk 
cylinder and the hollow tube behave quite differently. For the 
hollow tube the magnetic field is not expelled from the hole 
in the tube. The field remains as quantized flux with value 
quite close to that expected from a normal metal tube. The 
argument in Sec. II cannot explain this experiment. 

(5) As we stated in (4) the bulk superconductor can be re- 
placed by a superconducting solenoid whose leads are con- 
nected together The supercurrent flows in the solenoid and 
cancels the magnetic field inside. We can easily calculate the 
value of the current flowing. For a solenoid with winding 
number n per unit length, /o = B/tifio. Now what happens 
when the superconductor connecting the leads is replaced by 
a normal conductor This replacement can easily be done by 
raising the temperature locally. We do the operation at time 
f = 0. In this case, the current / decays in time according to 
the equation of the classical electrodynamics. 
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L— +RI^O, 

dt 



with the solution 
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(9) 



(10) 



where L is the self-inductance of the solenoid, and R is the 
resistance of the normal part. As the current decays, the mag- 
netic field penetrates inside the solenoid. At the same time 
Joule heating occurs in the normal part. The total amount of 
the Joule heat generated is easily calculated. It is same as the 
amount of the magnetic field energy that is introduced inside 
the solenoid. ' ' According to the argument in the previous sec- 
tion, the energy of the magnetic field is increased, since the 
solenoid is filled with the field. At the same time the same 
amount of the energy is dissipated as Joule heat. Where does 
this energy come from? How can the state realized as the 
minimum energy state change spontaneously to a higher en- 
ergy state? We cannot answer these questions, if we accept 
the argument in the previous section. 



IV. THE ORIGIN OF THE FLAW 

As we have seen in the previous section, the argument by 
Essen and Fiolhais brings consequences incompatible with ex- 
periments. Here we explain what is wrong in their argument. 
We need only knowledge of the electromagnetism to pinpoint 
their flaw. 

The essence of their argument is that the Meissner effect 
occurs to lower the magnetic field energy Eg in the volume of 
the conductor They believe that this energy can be reduced 



by flowing current at the surface of the conductor. However, 
this is wrong. As we stated in Sec. I, we need energy to cancel 
the magnetic field in the conductor. We calculate the energy 
needed to remove the magnetic field in several ways, and show 
that it is positive. 

We consider a situation in which a magnetic flux density B 
is created by an external magnet. The field is assumed to be 
almost uniform in the space in which we do experiments. We 
bring a solenoid whose length is I a, the cross section is S and 
winding number per unit length is n. We use this solenoid 
several times, so we name it as solenoid A. We place this 
solenoid A parallel to the external field B, and try to cancel 
the magnetic flux density inside the solenoid, the volume of 
which is y = IaS . It is noted that when the current flowing 
in the solenoid is zero, inner space of the solenoid contains 
magnetic energy 
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-B^y. 



(11) 



Let us assume for simplicity that the solenoid A is made of 
a perfect conductor, so the resistance is zero. The self induc- 
tance of the solenoid A is 



L^IJon^V. 



(12) 



In order to cancel the magnetic field inside, we need to flow 
current /o - B/fion. This current does not flow spontaneously. 
We need to apply electromotive force £(f) to the solenoid. The 
equation for the current is 



&(t)^L-I(t). 
at 

The work done by the electromotive force is 



(13) 
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(14) 



Here we used the condition that I(ti) - 0, and I(t2) - /o = 
B/tijj^). The amount of energy required is the same as that of 
the magnetic field stored in the volume V, but energy must 
be supplied externally. Namely, in order to remove magnetic 
field energy Eg from volume V, we need to input the same 
amount of energy to the solenoid. Where does sum of these 
energies, B^/fio = W + Eb, go? We will explain the destiny of 
this energy later. 

We can calculate the energy another way. A magnetic field 
exerts force on an electric current. Thus, the current in a 
solenoid is pushed by the magnetic field. This force is known 
as Maxwell stress, and the wall of the solenoid is pushed by 
pressure B^ IIjiq, the direction of which is from the high field 
side to low field side. Now we consider an infinitesimally 
thin superconducting solenoid placed parallel to the external 
magnetic field B. Since the self inductance of the solenoid is 
infinitesimally small, we can flow current to cancel the field 
inside without energy. Suppose that the solenoid is made of 



plastic material, and we can expand the radius while keep- 
ing the field inside to be zero. Since the wall of the solenoid 
is pushed from outside by the pressure P - B^/2jj.o, work is 
needed for this expansion. The work to realize a solenoid with 
inner volume V is PV = (B^/2yUo)V, the same work as before. 

We have shown that even though a finite amount of mag- 
netic energy is stored in space with density B^/2jj^), we need 
the same amount of energy to eliminate this energy. For the 
Meissner effect, this energy is supplied by the condensation 
energy of the superconducting state. We cannot extract the 
magnetic energy and utilize it without supplying extra energy. 
Essen and Fiolhais erroneously argued that the Meissner effect 
occurs to reduce this energy Eb without work W. As we have 
seen in Sec. Ill, their argument brings conclusions incompati- 
ble with experiments. 

Finally, we explain the energy balance. To eliminate ex- 
ternal magnetic field B from the volume V, we need work 
W = (B^/2yUo)y. This work and the energy of the magnetic 
field Eb must be absorbed somewhere. We can show that the 
energy is absorbed in the electromagnet that is creating the ex- 
ternal magnetic field B, or as potential energy of the magnet 
when the field B is due to a permanent magnet. For simplicity, 
let us consider the case in which the field B is created by an- 
other solenoid (which we label solenoid B) in which current 
Ib is flowing. 

The field at solenoid A should be proportional to Ib, so we 
write the relation as B = als- Since the magnetic flux that 
penetrates solenoid A is <1) = uIaSB - uviIaSIb, the mutual 
inductance between solenoid B and solenoid A is calculated to 
be M = O/Zb = ohIaS - anV. In the course of our increasing 
the current Ia in solenoid A to cancel the magnetic field B, 
electromotive force &b is induced in solenoid B>i^ 



diA dh 

dt dt 



(15) 



The work done by this electromotive force on the constant 
current Ib is 
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EiBlddt — anVlBh) , 



(16) 



where lA{t\) - and lAih) - h is used. Now we use 
the relation uIb - B and n/o - B/juq. Then we obtain 
Wb = (B^//V())y, which is just sum of the work we used to 
cancel the magnetic field in the solenoid W, Eq. (fT4l i. and the 
magnetic energy Eb, Eq. ( fTTT ). which existed in the volume. 
Our argument is consistent with classical electrodynamics. 

The situation in which the magnetic field inside the solenoid 
A is canceled by the current /q is a higher energy state, as we 
have emphasized several times. When we make the solenoid 
A to have finite resistance R, the work W spent to realize the 
situation is returned from the source of the magnetic field, i.e. 
from the power source of the solenoid B, as Joule heat gener- 
ated in solenoid A. Actually, it is calculated as follows, 

Rlitfdt^ I 7;/2exp(-2-fj 
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(17) 



where I(t) = Io[-{R/L)tl Eq. ^, L = fion^V, and /q = 
B/n/io are used. 

Since the state where the magnetic field is expelled has 
higher magnetic field energy, the Meissner effect is possible 
only when the condensation energy to the superconducting 
state is larger than the magnetic field energy. This condi- 
tion determines the critical field H^. For ordinary metals the 
magnetic susceptibility is quite small, so the free energy of 
the normal metal, F„ is almost independent of the magnetic 
field applied. If the free energy of the superconducting state 
in the absence of the magnetic field is F^, the critical field 
He = Bc/i-iQ is determined by the condition^^iii 



F, + -yU0//2y = Fn . 



(18) 



In this equation V is the volume of the superconductor, and 
F„ is the free energy when the superconductivity is destroyed 
by the magnetic field. This relation has been experimentally 
confirmed. 



V. CONCLUSIONS 

As is commonly understdoood, the Meissner effect requires 
quantum mechanics for its explanation. We have reviewed 
established evidence for this in the introduction. This fact 
was challenged by Essen and Fiolhais in a recent publication. 
We have explained their argument in Sec. II, and have seen 
in Sec. Ill that if we accept their argument, unphysical phe- 
nomena are anticipated. We clarified what is wrong in their 
argument. 
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